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1. INTRODUCTION 

The study of statistical physics models in two dimensions (d = 2) at their critical point 
is in general a significantly hard problem (not to mention the d = 3 case). In the eighties, 
three physicists, Knizhnik, Polyakov and Zamolodchikov (KPZ) came up in [KPZ88J 
with a novel and far-reaching approach in order to understand the critical behavior of 
these models. Among these, one finds for example random walks, percolation as well as 
the Ising model. The main underlying idea of their approach is to study these models 
along a two-step procedure as follows: 

- First of all, instead of considering the model on some regular lattice of the plane 
(such as Z 2 for example), one defines it instead on a well-chosen "random planar 
lattice". Doing so corresponds to studying the model in its quantum gravity form. 
In the case of percolation, the appropriate choice of random lattice matches with 
the so-called planar maps which are currently the subject of an intense activity 
(see for example [LGMllb] ). 

- Then it remains to get back to the actual Euclidean setup. This is done thanks to 
the celebrated KPZ formula which gives a very precise correspondence between 
the geometric properties of models in their quantum gravity formulation and their 
analogs in the Euclidean case. 

It is fair to say that the nature and the origin of such a powerful correspondence 
remained rather mysterious for a long time. In fact, the KPZ formula is still not 
rigorously established and remains a conjectural correspondence. The purpose of this 
survey is to explain how the recent work of Duplantier and Sheffield enables to explain 
some of the mystery hidden behind this KPZ formula. To summarize their contribution 
in one sentence, their work implies a beautiful interpretation of the KPZ correspondence 
through a uniformization of the random lattice, seen as a Riemann surface. 
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To fully appreciate the results by Duplantier-Sheffield, we will need to introduce 
beforehand several related concepts and objects. More precisely, the rest of this in- 
troduction is divided as follows: first we give a short and informal discussion about 
quantum gravity, then we introduce two universality classes of random lattices. Then 
we will come to the KPZ formula through a specific example (boundary of Random 
Walks hulls), and finally after stating the main Theorem by Duplantier-Sheffield, we 
will state a beautiful conjecture they made. 

1.1. A first glance into quantum gravity 

Quantum gravity is intimately concerned with the following naive question: 

QUESTION 1.1. - - How does a "uniformly distributed" random metric on the sphere S 2 
typically look ? 

What is naive in this question is the fact that one would first need to specify what 
we mean by "a uniform probability measure" on the space of metrics on § 2 . It turns out 
that defining a natural model of random metric on the sphere § 2 already is a difficult 
and interesting problem. To illustrate this, let us ask a similar naive question in a 
one-dimensional setting: 

QUESTION 1.2. - - For any a,b G M d ; how does a "uniformly distributed" path 7 : 
[0, 1] lR d going from a to b typically look ? 

This naive question was of crucial importance at the time Feynman developed the 
so-called path integral formulation of quantum mechanics. 

Already in this case, defining properly a "uniform measure" on paths was not an easy 
task. Yet, it had been mathematically settled prior to Feynman's work and corresponds 
to the well-known Brownian motion. 

In some sense, the purpose of quantum gravity is to extend Feynman path integrals to 
Feynman integrals over surfaces. Physicists are particularly interested in such an exten- 
sion, since this would provide a powerful tool to deal with the quantization of gravitation 



field theory, a notoriously hard problem. W With this background in mind, the problem 
of defining a proper mathematical object for a "uniformly chosen random metric on § 2 " 
thus corresponds to defining a two-dimensional analog of Brownian motion, i.e. a kind 
of Brownian surface. 

Even though physicists are primarily interested in the above continuum formula- 



tion of Question 1.1, a natural and very fruitful approach is to study an appropriate 



discretization of it and then to pass to the limit. This brings us to the next subsection. 



1. This approach towards the quantization of gravitation is called loop quantum gravity. 
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1.2. Discretization of Question 1.1 and planar maps 



In the one-dimensional setting, if one asks Question 1.2, it is not straightforward to 



come up right away with Brownian motion. But, if instead we start by discretizing 



Question 1.2, say by allowing random 1/n-steps, then we end up with the model of 
random walks. Brownian motion is then obtained as the scaling limit (as n — > oo) of 
these rescaled random walks. 



It is thus tempting to apply the same strategy to Question |1.1[ namely to find an 
appropriate discretization. Let us explain below a possible discretization which was 
used initially in the physics literature and was studied extensively recently among the 
mathematical community. See for example [LG Mllb] and references therein. We will 
see in subsection |1.3 that there are other ways to discretize Question LI which lead 



to different universality classespH but the discretization below is in some sense the 



simplest and most natural one regarding the statement of Question |1.1 

The idea of the discretization we wish to introduce is to consider discrete graphs, 
with say n faces, which have the topology of a sphere S 2 and for which the metric p n 
will correspond (up to a rescaling factor) to the graph distance, i.e. p n := n~ a d gr for 
some exponent a > 0. The exponent a will need to be well chosen as in the case of 
Random Walks, where space needs to be rescaled by y/n in order to obtain a limit. 
If we define our discretization in such a way that for each n > 1, there are finitely 



many such graphs, we can pick one uniformly at random (in the spirit of Question 1.1 ) 
which thus gives us a random metric space (M n ,p n ). We can then ask the question of 
the scaling limit of these random variables (M n ,p n ) as n — > oo in the space (K, daii) 
of all isometry classes of compact metric spaces, endowed with the Gromov-Hausdorff 



distance oIgh ^ The advantage of this setup is that (IK, den) is a complete, separable, 
metric space (a Polish space) and is thus suitable to the analysis of convergences in law 
and so on. Note here, that even if one could prove that (M n , p n ) converges to a limiting 
random object (Moo,Ax>)) it is not clear a priori that the topology is preserved at the 
scaling limit or in other words, it needs to be proved whether (Moo, Poo) a.s. has the 
topology of a sphere or not. If all these steps can be carried on, then this would give 
us a good candidate (Moo, Poo) for the random object used in Question |1.1| 

Let us now introduce one specific discretization. 

Definition 1.3 (planar map, following |LGllj ). — A planar map M is a proper em- 
bedding of a finite and connected graph into the two-dimensional sphere S 2 , which is 
viewed up to orientation preserving homeomorphisms of S 2 (i.e. up to "deformations"). 
Loops and multiple edges are allowed. The faces of M are identified with the connected 



2. Similarly as Random Walks with non-L 2 steps converge to other Levy processes than Brownian 
Motion 

3. Informally, if (Ei, di) and (E 2 , d 2 ) are two compact metric spaces, then d,Qfj{Ei 1 E 2 ) is computed 
as follows: we embed E\ and E 2 isometrically into some larger metric space (E 7 d) and we compute 
using the common distance d the distance in the Hausdorff sense between the two embeddings. Then 
we take the infimum over the possible such embeddings. Sec [LGMllb . 
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components of S 2 \ M and the degree of a face f is defined as the number of edges 
incident to f , with the additional rule that if both sides of an edge belong to the same 
face, this edge is counted twice. 

Finally, for combinatorial reasons, it is often convenient to consider rooted planar 
maps, meaning that one particular oriented edge e is distinguished. The origin of that 
root edge it is called the root vertex 0. See figure [l] for an instance of a planar map 
where all faces happen to be squares. 



FIGURE 1. This is a planar map of the sphere S 2 with exactly 17 squares (this 
includes the exterior square which is also in the sphere). 

Definition 1.4 (p-angulations of the sphere). — For any integer p > 3, let be the 
set of all rooted planar maps with n faces, where each face has degree p. The elements 
of are called rooted planar p-angulations. (p = 3 corresponds to triangulations 
and p = 4 to quadrangulations). 

Since planar maps are defined up to deformations, there are finitely many rooted 
planar maps in M£, for any n > 1. For example, one has 



The appearance of the n-th Catalan number here is explained through the celebrated 
Cori-Vauquelin-Schaeffer bijection, which gives a one-to-one correspondence be- 
tween labelled plane trees and rooted quadrangulations of the sphere (see for example 





[LG MTTb] ). 
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For a fixed integer p > 3, we will denote by m„ G a sample of a planar map 
uniformly distributed over M^. Since in Question we were interested in a uni- 
formly distributed random metric on S 2 , it makes sense to consider the random variable 
(m n , n~ a d gr ) seen as a random point in the above space (K, (Igh)- 

The study of these random planar maps has now a long history (see for example 
jLGMllbl ILGP081 ILG071 IMieOSj ) and has culminated in the following breakthrough 
result proved independently by Le Gall QLGllJ) and Miermont ( jMiellj 



(4) 



THEOREM 1.5 (Uniqueness of the scaling limit). - - There exists a (unique) random 
compact metric space (m^, D*) with values in K such that for any p G 3 U (2N + 4), 
there is a positive constant X p such that as n — )■ oo, 

( mn ' ~i74 d Sr) ( m °°' D *) ' 

in the Gromov-Hausdorff sense. This random compact metric space is called the Brow- 
nian map. 



Furthermore the following property holds. (It was first established in |LGP08] . See 
also |Mie08] for a different proof.) 



THEOREM 1.6 (Sphericity of the Brownian map). — Almost surely, the random met- 
ric space (moo,!)*) is homeomorphic to the sphere § 2 . 



This last theorem thus ensures that in some sense the Brownian map (m^, D*) is a 
good candidate for the "uniform" random metric on § 2 considered in Question |l.l[ The 
only drawback being that the metric space (moo,/}*) which is a.s. homeomorphic to 



(S> 2 , || • ||r3) is not provided with a "canonical" embedding into the sphere^ We will 
come back to this question of embedding later. 



1.3. Quantum gravity coupled with "matter" 

Coming back to our earlier motivation, i.e. the study of statistical physics models on 
regular lattices, we will now introduce models of random lattices which are naturally 
associated with statistical physics models. The first one will be associated to the Ising 
model and will fall into a different quantum gravity universality class than the above 
planar maps. 



4. To be more precise, Miermont's proof is restricted to the case p = 4 but gives slightly stronger 
estimates on the structure of the geodesies while Le Gall's proof is more general in the sense that it 
proves the result also for p = 3 and all even p greater than 4. 

5. This would be the case for example if (moo, D*) happened to be a nice Riemann surface, but it 
is not. 



1052-06 



1.3.1. Quantum gravity coupled with Ising model. - - The standard and simplest way 
in statistical physics to model ferromagnetic matter (like a piece of iron for instance) 
is through the so-called Ising model. It is defined on a graph G = (V, E) which is 
supposed to represent the metallic structure of our ferromagnet. For a flat piece of 
iron, one might choose Z 2 for example. On a finite graph G = (V,E), it is defined as a 
probability measure on spin configurations a G { — 1, where if x G V, a x — +1 
means that the spin of the atom at site x is oriented in the upward direction. The 
probability measure P = Pj is such that each spin configuration a has a probability 
proportional to 

Pj [a] oc exp ^ J cr x a y j , 

e=(x,y) 

where the parameter J represents the strength of the electromagnetic interaction be- 
tween atoms . Note that the higher J is, the greater the tendency is for the spins 
to align in the same direction. In order to express the probability measure P more 
explicitly, it is natural to introduce the partition function of the graph G 

Zj = Zj(G):= ^2 exp^J ^2 a * a 

so that 



y 

ae{-l,l} v e=(x,y) 



v , y 

e=(x,y) 

The Ising model has been studied extensively in the physics and mathematical com- 
munity over the last 70 years. See for example |SmilO| and references therein for a 
recent breakthrough paper on the subject. 

Now the idea of studying such a model from a quantum gravity perspective is to 
sample a spin configuration a together with a random base-graph G. This is what one 
calls coupling quantum gravity with the Ising model (i.e. we are looking for a probability 
measure on pairs {G = (V,E),a G { — 1, 1}^} )• One way to proceed would be to first 
sample the random lattice G, using the above planar maps m n uniformly chosen among 
and then to sample an Ising configuration a on the graph m n according to the above 
model P = Pj,m n . Doing so would not correspond so much to coupling quantum gravity 
with the Ising model, since the marginal of (m n , a) on the graph- component would just 
be the above model of planar graphs. So we will use a more intricate way of coupling 
the graph with its spin-configuration. 

If we restrict ourselves to the case of planar graphs, the standard quantum grav- 
ity/Ising coupling works as follows: fix p > 3 and some large n > 1. For integrability 
reasons, the spins will be indexed by the faces of the planar maps m G instead 
of the vertices. (In other words, we will consider our Ising model on the dual graphs 



6. Often, instead of the parameter J, one uses the inverse temperature j3 := -^Kp and the Ising 
model is defined so that P/3[c] oc e l3 ^"=^-v) a " av ; but mainly for notational reasons, we will not use 
this point of view here. 
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m* which are no longer in M^). With this slight change, a natural idea is to define a 
measure P = Pj >n on couplings (m n ,er) £ x { — 1, l} n (recall that planar maps in 
have exactly n faces) so that 

F J>n [m, a] oc exp (j ^ a f a f') > 

/-/' 

where the sum is over pairs of adjacent faces in m. See Figure [2} 




FIGURE 2. This is an Ising spin configuration on the same quadrangulation as 
in Figure [T] (were we removed the inside edges for convenience) . 

Note that, unless J = 0, the marginal on m (the graph-component) is no longer 
uniformly chosen among M^, but follows a law which highly depends on the magnetic 
strength J. In probabilistic terms, the marginal on m n corresponds to weighting m n 
by its Ising partition function Zj(m n ). As in the case of a fixed graph G, it is natural 
to define the following annealed partition function 

Zn{J) ■= Yl eXP ( J ^2 a f a f') = Yl Z A m n) 

meM£,cre{-l,l} m m6M£ 

Now, recall the KPZ formula is supposed to relate a critical Euclidean Ising model 
(say on Z 2 ) with a critical Ising model defined together with a random planar map 
m G M in the above fashion. Therefore, we need a way to detect a phase transition 
in each model, Euclidean and quantum gravity. Among the many ways to "detect" phase 
transitions, one is to notice a failure of analyticity for certain thermodynamic quantities. 
In our case, partition functions enable us to detect a phase transition. More precisely, 
for the Euclidean Ising model on Z 2 , for each n > 1, let Z n (J) be the above partition 
function for the graph G n = [— n, n} 2 (1 1? . Then, from the work of Onsager |Ons44j . it 
can be shown that as n goes to infinity, the functions J \— > log Z n (J) converge to a 
limiting function f(J), called the free energy which is analytic for all values of J > 
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except at J = J° ncl = | log(l + \/2). On the quantum gravity side, the natural quantity 
to look at is the so called grand canonical ensemble partition function defined as 

n>l meM%,<re{-l,l} n 



n>l 



The reason why such an annealed partition function is considered is because it was 
computed exactly by Kazakov ( |Kaz86j ) in the case of quadrangulations (p = 4). His 
computation relied on a deep relationship with some (random) matrix models. This 
type of exact computation occurs for what one calls exactly solvable models and is the 
reason why we changed slightly the definition of Ising model from vertices to faces. The 



exact computation by Kazakov enables to detect a bi-critical point (/3 C , J c ) ^ around 
which the analyticity of Z(/3, J) is broken. The analysis of Z around its bi-critical point 
gives detailed information on the critical behavior of the system in its quantum gravity 
formulation. See |LGMlla] and Appendix B in |Dup04| for thorough discussions of 
this. 

As mentioned above, the scaling limit as n — > oo of these planar maps m n G 
weighted by Zj c= i n 2(m n ) will fall in a different universality class than the "uniform" (or 
unweighted) planar maps. We will come back to this later. 

1.3.2. Quantum gravity coupled with random walks, self avoiding walks or percolation. 

- If one wants to study statistical models such as random walk, percolation or self- 
avoiding walks (SAW), it turns out that all of them are naturally coupled with the 



planar maps introduced in subsection |1.2[ This regime corresponds to what one calls 
the pure gravity regime. It means in some sense that the random geometry of the 
planar map is insensitive to the model it is coupled to (as opposed to what happens 
with Ising model). We will describe this pure gravity case in more details in the next 
subsection through the example of random walk. 



1.4. The KPZ formula in the special case of random walk and Brownian 
motion 

Let us start by the Euclidean case. Consider a random walk X t on the domain 
A A? := [-N, N] 2 starting at the origin until it reaches the boundary 8An- One can look 
at several subsets of interest about this random walk among which: 

(i) The range IZn of X t , i.e. IZn '■= {X t }. 

(ii) The set of cut-points Cn of X t , i.e. the set of points x G An so that removing x 
disconnects the range into 2 disjoint components. 

(iii) The set of frontier points J^n, i-e. all the points on the range that are connected 
to 8An in the complement An\7Zn- 



7. When p = 4, one has J c = J® G = In 2 
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See figure [3j One way to measure the typical size of a random subset of A^r is via the 
following notion of scaling exponent. If K = (Kn)n>i is a certain sequence of random 
subsets of A^v, then its Euclidean scaling exponent is defined as 

logE[|fGv|/iV 2 ] 



(1) 



x 



x(K) :-- 



lim 



logl/N 2 

Note that we are being informal here since it might be that this limit does not exist. 
To be more rigorous, one should define instead upper and (lower) scaling exponents 
by using lim sup (and liminf) instead. Nevertheless, to simplify the exposition we will 
neglect this issue in the remaining of this text and will assume that the random sets 
we will consider are such that these limits always exist (proving such a convergence 
can be very hard for some models and remains in many cases an open problem, for 
example scaling exponents for critical percolation on Z 2 ). Note also that the greater 
the exponent x = x(K) G [0, 1] is, the "smaller" the random sets are since having 
scaling exponent x means that asymptotically, is of size ~ N 2 ~ 2x . 



A 



N 



m 



5 



X 




Figure 3. 

The scaling exponents for the above subsets of random walk paths have been 
the focus of an intense activity over the last 20 years or so. In fact the road which 
culminated in their evaluation is quite an interesting story. Their values were first 
conjectured by Duplantier and Kwon in [DKH88J based on conformal invariance and 
numerical simulations. Somewhat surprisingly, the powerful machinery of conformal 
field theory was unable at giving theoretical predictions for these exponents. In 1998, 
Bertrand Duplantier obtained in |Dup98| (on a non-rigorous basis) these exponents 
with a completely different approach: his idea was to study a random walk path on 
a random planar map m ra and to compute the scaling exponents, called quantum 
scaling exponents for the frontier and cut points of the random walk on m n . Once this 
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was achieved, he was able to recover his previous conjecture via the KPZ formula. His 
work was partly based on the seminal work [LW99J on Brownian intersection exponents. 
Finally, the story ends with the mathematical derivation of these exponents in |LSW01a] 
LSWOlbJ which used yet another approach: the Schramm-Loewner-Evolution SLE K 
processes (with k = 6 here). 

Here is more precisely how it goes. Imagine that we have at our disposal a natural 
model of planar maps m n with boundary (i.e. with the topology of a disk), instead of 
planar maps m n with the topology of a sphere. This can be done in such a way that 
it gives a different, yet interesting model of planar maps: see |Betllj where Bettinelli 
considers random quadrangulations m n with n faces and which are such that dm n is of 
order yfn. Let thus m n be sampled uniformly among rooted quadrangulations with n 
faces and boundary of length y/n. Let W m be a simple random walk on this graph 
m n starting at the root of m n until it reaches the boundary dm n . See figure [3j One 
can consider the exact analogs of the above sets, i.e: 

(i) The range TZ n of W m , i.e. Tt n : = {W m }. 

(ii) The cut-points C n of W m , i.e. the set of points x G m n so that removing x 
disconnects the range into 2 disjoint components. 

(iii) The frontier points J-" n , i.e. all the points on the range that can be connected to 
dm n via the complement m n \ TZ n . 

Similarly to the Euclidean case, if K = (K n ) n >i is a sequence of random subsets of 
m n , then its quantum scaling exponent is defined as: 



(2 ) A = A(A") := U m gM^j = Um gliggWH , 

n^oo logl/|m n | n->oo logl/Tl 

where, once again, we implicitly assume here that the limits exist. 

Even though it was not exactly in the same setup, Duplantier determined in |Dup98 
the values of these quantum scaling exponents. He found A (J 7 ) = 1/2 and A(C) = 3/4. 
Now for this universality class (i.e. the pure gravity case of "uniform" planar maps 
m n ), the KPZ formula reads as follows: the scaling exponents x and A of a certain 
subset of a random walk path in the Euclidean and quantum cases are related by the 
following quadratic expression: 

(3) s = ?A 2 + ±A. 

One thus finds x(T) = 1/3 and x(C) = 5/8. Note that for the range 71, it is well 
known that an Euclidean random walk X t visits about N 2 / logiV sites before touching 
<9An, this means that x(7V) = which is a fixed point of the above KPZ formula. This 
thus suggests that A {TV) = as well, which in turn means that the random walk W m 
should visit iV^m™)! 1-0 * 1 ) = n l ~°^ sites of the planar map m n . This is a non-trivial 
fact in its own, which illustrates that the KPZ formula (if true) can be powerful in both 
directions. 
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1.5. The KPZ formula in general 

The planar maps m n uniformly chosen among some for some p > 3 are the natural 
universality class for studying random walk on the quantum gravity side. Nevertheless, 
one might also consider an independent random walk on a planar map m n sampled 
together with an Ising model as we explained above. If one would do so, this would 
affect the typical size of the sets C n and T n and one would find different values for 
the quantum scaling exponents A ising (J-') and A ising (C). This reflects the fact that the 
two different procedures we gave to sample planar maps do not fall asymptotically in 
the same universality class. For this Ising universality class, the KPZ relation reads as 
follows: 

(4) x — ^^ising + ^ Rising • 



Such a correspondence can be very useful to relate critical exponents for the Ising 
model in its quantum gravity form with the analogous Euclidean critical exponents. For 
example, one might hope that if K n denotes the largest + cluster in the random map m n 
weighted according to Zj(m n ), then its quantum scaling exponent Ai S ; ng should relate 
to its Euclidean analog through the above formula. Yet, one has to be careful with 
such statements since, especially away from the pure gravity case (i.e. uniform planar 
maps), it is in general a subtle affair to know which critical exponents are allowed to 



go through the KPZ formula or not. W In the case of the Ising model, quantum critical 
exponents are usually derived from the analytic study of Z(/3, J) near its bi-critical 
point and are then transfered to the Euclidean setting using Q. See |Dup04| . 

One can now state the general KPZ formula: consider some statistical physics model 
defined together with a random planar lattice (in the suitable universality class). Then, 
the quantum and Euclidean scaling exponents A and x which describe the critical 
properties of the model, are such that the following quadratic KPZ formula holds 

(5) x = ^A 2 + (l-^)A, 

where the parameter 7 G M+ determines in which universality class the model is. We 
thus find that pure gravity corresponds to 7 = \/8/3 while Ising model corresponds 
to 7 = \/3. More generally, it follows from this correspondence that critical two- 
dimensional models form in some sense a one-dimensional family space. This is consis- 
tent with the recent theory of SLE K processes which are aimed at describing interfaces 
of critical two-dimensional systems and also form a one-dimensional family of processes 
indexed by k > 0. In fact, the two parameters are related one to another by 7 = yfn. 



8. The subtlety lies partly in the fact that one needs the set K to be independent of the field e lh in 
the main Theorem 1.10 while here the largest + cluster would be highly correlated with the field e lh . 
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1.6. Why go through quantum gravity ? 

The KPZ formula becomes particularly useful (at least on a non-rigorous level) when 
the critical exponents of a particular model happen to be much easier to compute, say, 
in its quantum gravity form than in its Euclidean one (or vice- versa). This is exactly 
what happens in the case of random walk ( |Dup98| ). 

1.6.1. Quantum gravity coupled with random walk. - - . Imagine we want to study the 
exterior frontier Tt of a random walk X t in Z 2 as t increases (J- t is defined here as the 
set of points in lZ t which are connected to infinity in Z 2 \ 1Z t ). For each T > 1, the 
evolution of Tt for t > T will not depend of what the random walk did inside the hull 
T>t, defined as the complement of the unique infinite connected component of Z 2 \ TZt- 
Yet the evolution of Tt,t > T will depend of the complicated law of the boundary 
8T>t = Tt C Z 2 which as T — > oo looks more and more like a fractal set. On the 
quantum gravity side, such geometric processes as Tt behave in some sense in a much 
more Markovian way. 

To illustrate this, let us briefly describe a striking recent work by Benjamini and 
Curien [BC12J. In their work, they rely on an infinite version of the above rooted 
quadrangulations of the sphere m n G M^. Qoo is called a rooted UIPQ (for Uniform 
Infinite Planar Quadrangulations) and can be seen as a local limit of the quadrangu- 
lations m n as n — > oo. They consider a random walk W m starting at the root of Qoo 
and they explore the planar map Qoo along the random walk path W m . The crucial 
observation which goes back to |Ang03| is that if at time m, T m denotes the exterior 



boundary of the domain which was explored by the random walk so far^\ then the law 
of the infinite connected component Qoo \ T> m depends only on the length of dT> m = T m 
and not at all on its shape. This simplifies things tremendously. Yet, it would remain to 
control how \T m \ and [D m \ both increase as a function of m. This part turns out to be 
difficult but Benjamini and Curien are able to obtain in |BC12j that as m — > oo, \T m \ 
behaves like iVm] 1 ^ 2 . This provides a rigorous derivation of the identity A (J 7 ) = 1/2. 
(They also obtain another quantum scaling exponent for the so-called frontier points). 

These nice spatial Markovian properties were also crucial in the seminal work on 
the quantum scaling exponents of random walk by Duplantier (|Dup98|), where various 
quantum scaling exponents for the random walk were determined through an asymptotic 
analysis of the partition function of random walks coupled with planar maps. 

1.6.2. Quantum gravity coupled with percolation. - - Critical percolation has also been 
successfully analyzed from the quantum gravity perspective. Various quantum criti- 
cal exponents for percolation were determined in [SD87J and were translated to the 
Euclidean setting via the KPZ formula. This is probably the first explicit use in the 
physics literature of the KPZ formula. In the mathematics literature, there is a work 



9. some care is needed here, since one has to prove that under the law of Qoo, there is a unique 
infinite connected component in Qoo \ T^mi but we will not enter in more details here. See |BC12j . 
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in progress by Angel and Curien [ACJ, which gives among other things a rigorous de- 
termination of the critical points for bond percolation on the UIPT and UIPQ (the 
first one is the analog of the above with triangles instead). Using the Markovian 
structure we briefly sketched above, they are able to compute some of the quantum 
scaling exponents of these critical percolation models, which is a very interesting step 
towards the understanding of statistical physics models on random lattices. 

1.6.3. Quantum gravity coupled with Ising model. - - The reason why it is interesting to 
study the Ising model in its quantum gravity form lies in the fact that (as was mentioned 
above), Kazakov was able to compute exactly the annealed partition function 

n>l meM* 

As hinted previously, various quantum critical exponents for the Ising model can be ex- 
tracted from the behavior of 2(8, J) near (8 C , J c ) and can be translated to the Euclidean 
setting via the KPZ formula ^ with 7 = y/3. 

It is now time to introduce the main result we wish to explain in this survey. 



1.7. The main result by Duplantier-ShefReld and a conjecture on the em- 
bedding of planar maps 

What remains mysterious so far in the KPZ correspondence we just described is the 
fact that it relates scaling exponents of sets which do not "live" on the same space. 
For example it relates the scaling exponent of cut-points of a random walk living in A^v 
with the (quantum)-scaling exponent of cut-points of a random walk living on a random 
planar map m„. To overcome this, Duplantier and Sheffield discovered a setup in which 
the random set we are interested in, say K, lives on some single space (say [0, l] 2 or S> 2 ), 
but its size, given in terms of its scaling exponent, can be measured in two different 
ways: an "Euclidean" way which will output some Euclidean scaling exponent x = x(K) 
and a "quantum" way which will output a possibly different quantum scaling exponent 
A = A(K). Their setup is built so that x(K) and A(K) will satisfy to the general KPZ 
relation (|5]). (More precisely for each value of 7 > 0, they found an adequate setup). 

Before going into the details of their setup, one sees here that we are getting closer 
to our initial naive Question |1.1[ Indeed, if one had at our disposal a "uniform" random 
metric p on S> 2 , one could evaluate the size of K either using the Euclidean metric 
|| • || (which would give us a scaling exponent x) or using the random metric p (which 
would give us an exponent A). Now from the preceding discussions, x and A should 
satisfy to (5) with 7 = ^/8/3. Unfortunately we did not quite succeed yet in answering 
Question [LT] Indeed, Theorems L5 and [L6] provide us with a "uniform" probability dis- 
tribution on compact metric spaces which a.s. have the topology of a sphere, but there 
is (at least for the moment) no canonical way to embed the limiting space (moo,/}*) 
into the sphere. In fact, because of this embedding issue, a proper answer to Question 



1.1 remains an important open problem. 
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Yet, before passing to the limit n — > oo, there are several ways to "naturally" embed 
the planar maps m n into the sphere. We present two of them. 

1.7.1. Embedding of planar maps seen as Riemann Surfaces. - - The main idea here 
is that one can view each planar map m n G as a Riemann surface. For this, view 
each face as a polygon (here a square) on which we give the obvious flat conformal 
structure given formally by z i— > z. By the Schwarz reflexion principle, one can easily 
glue together the conformal structures of two adjacent faces along their edge. Some 
care is needed around each vertex x G m n since the angle around x might not be 2n. 

4 

Yet, one can use a local chart around each such x of the form z t- >• zk where k is 
the degree of vertex x (this corresponds to the fact that conic singularities of complex 
manifolds are removable). Altogether, this gives us a complex manifold of dimension 
one endowed with a finite atlas indexed by the set of edges and vertices of m n . By 
the Riemann uniformization theorem, since m n ~ S> 2 , m ra equipped with the above 
conformal structure can be mapped conformally to § 2 and the embedding is unique 
up to Mobius maps of the sphere § 2 . We thus found a natural way (up to Mobius 
transformations) to embed any planar maps m„ into § 2 . The same idea enables us to 
embed in a conformal way planar maps m n with the topology of a disc (with dm n ^ 0) 
into the disc D or into [0, l] 2 . 

1.7.2. Embedding of planar maps via circle packings. - - If one considers planar maps 
m n G (i.e. triangulations of the sphere), then by Kobe Theorem, there is a unique 
(again up to Mobius transformations) circle packing in § 2 whose connectivity graph 
corresponds to m n . See figure |4j 

This circle-packing embedding is different from the one given by the conformal struc- 
ture, yet one might conjecture that if m n is sampled uniformly from M^, then the two 
embeddings should look almost alike with large probability. 

1.7.3. Scaling limit of planar maps embeddings. - - Let m n C § 2 be a "natural" embed- 
ding of a random planar map m n G into S 2 (m n may be sampled either uniformly 



as in subsection 1.2 or weighted according to some statistical physics model as in sub- 
section 1.3). To each embedding m n , the renormalized graph distance n~ a d gr on the 
vertices of m„ C § 2 can be easily extended to a metric p n on the whole sphere § 2 . 



If m n is sampled uniformly in M£, then we know from Theorem 1.5 that one should 
choose a = 1/4. If one could prove in this case that as n — > oo, the (random) metric p n 
on § 2 would converge in law, then it would give a "good" answer to Question |l.l| and 
it would enable us to build a setup for a concrete interpretation of the KPZ formula 
when 7 = a/8/3. 

In some sense the idea of Duplantier and Sheffield is to focus on measures instead 
of distances. If m n C § 2 is a natural embedding of a planar map m n , then it is natural 
to consider the pushforward in S 2 of the Lebesgue measure on m n , renormalized so 
that m n has unit area (i.e. all faces of m n have Lebesgue measure exactly n^ 1 ). Let 
us denote by \x n this pushforward measure on S 2 . fi n is thus a random measure on the 
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FIGURE 4. This picture is a simulation by Maxim Krikun which represents a 
circle packing of a uniform triangulation on the disc. 

sphere with /i n (§> 2 ) = 1. As one can see from Figure [4j one expects that //„ should look 
quite singular with respect to Lebesgue measure on § 2 (or D if one considers maps with 
the topology of a disc). Let /„ be the Radon-Nikodym derivative of fi n with respect to 
the Euclidean Lebesgue measure, then we expect f n to become more and more "rough" 
as n — > oo. This brings us to the question. 

Question 1.7. — 

(i) If m n is sampled uniformly in is it the case that fi n converges as n — )■ oo to 
a random measure on S 2 ? 

(ii) What if m n G is weighted by Z Jc (m n ) ? 
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Solving this question is a major open problem in the area, but Duplantier and 
Sheffield made a decisive step in this direction: they managed to identify an ex- 
plicit candidate for the scaling limit of /i n , for which the KPZ equation holds (in the 
sense of measures). This candidate as we will see in more details below is given by the 
exponential of a Gaussian Free Field. 

1.7.4- Setup and statement of the main Theorem. - - The common base space will be 
[0, l] 2 and we will consider some deterministic subset K C [0, l] 2 . Now, we consider 
two measures on this space: 

(i) The Lebesgue measure C on [0, l] 2 

(ii) A random measure /x = /i 7 which can be formally written as = e 7h , where h 
is an instance of a Gaussian Free Field. 

The Gaussian Free field (GFF for short) is a certain Gaussian process which will be 
defined in section [2] It is "too rough" to live in the space of functions on [0, l] 2 (say 
L 2 ([0, l] 2 ) ) and has to be viewed instead as a random Schwartz distribution. See Figure 
[5] for a representation of how a regularized GFF on D looks like. 




Figure 5. An instance of a Gaussian Free Field h in the disc D with Dirichlet 
boundary conditions. Picture by Nam-Gyu Kang. 

Since h is not a proper function, e lh is not well-defined a priori. This will be the 
content of section [3] to give a meaning to such measures. See figure [6] for an illustration 
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of these measures when 7 G {3/2,5}. These random measures are supposed to model 
the effect of quantum gravity on our base-space [0, l] 2 . They are called Liouville 
measures. 
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FIGURE 6. In each picture, all the plain dyadic squares have about the same 
quantum area 5. The picture on the left is for 7 = 3/2 while the picture on the 
right which looks much more singular is for 7 = 5. 



Let us now define Euclidean and quantum scaling exponents in this setting. For this, 
following the notations of [DSllJ, we will need the following notion of Euclidean and 
quantum balls. 

Definition 1.8 (Euclidean and quantum balls). — For any point z G [0, l] 2 and any 
e, 5 > 0, define 

(i) B e (z) to be the Euclidean ball of radius e around z. 

(ii) B 6 (z) to be the ball B T {x) with r := sup{r > 0, n(B T (z)) < 5}. B s (z) will be 
called the quantum ball around z of quantum area 5. 



(10) 



Definition 1.9 (scaling exponents). — Let K C [0, l] 2 be fixed. 
1. The Euclidean scaling exponent x = x(K) is defined as 

\og¥[B e (z) n K ^ 0] 



(6) 



x 



x(K) :-- 



lim 

€->0 



loge 2 



(7) 



where the point z is sampled uniformly on [0, l] 2 . 
2. The quantum scaling exponent A = A(K) is defined as 

\ogE[fi[B s (z)nK 7^0]] 



A = A{K) :-- 



lim 

<5->0 



log 5 



10. Note that since we do not have at our disposal any "quantum distance" on S 2 , our definition of 
quantum balls still relies somewhat on the Euclidean one. Nevertheless, it is believed that this slight 
Euclidean use should "average out" as S —> 0. 
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where z is "sampled" according to the a.s. finite measure a = e lh and E averages 
over the random measure //. 

As we did before, we implicitly assume here that the limits exist. 

Note that these definitions are the exact analogs of equations ^ and ([7]). Indeed if 
one believes in fi n — y \i and if one approximates K be a subset K n C m n , then with 
5 = 1/n one has 

a ^ iogE[^[g^n^0]] 

log 5 

log E [// n [ the face containing z ~ /z n in m n is contained in K n \ ] 

r~^j 

log 1/n 

logE[|Ar ra |/n] 
log 1/n 

and similarly for x = x(K). We are now in position to state the main theorem of 
Duplantier and Sheffield. 

THEOREM 1.10 (Duplantier and Sheffield [DSllJ). - - Consider the Liouville measure 
/x 7 = e lh on the unit square [0, l] 2 and let K be a (deterministic) subset of [0, l} 2 , such 
that the limits in |6p and |?p exist. Then if 7 G [0,2), the quantum scaling exponent 
A = A(K) for the Liouville measure /i 7 almost surely satisfies to the KPZ formula: 

(8) x = ^A 2 + (l-^)A. 

Remark 1.11. — Note that in the statement of the theorem, by Fubini's theorem, K 
may also be a random subset of [0, l] 2 but in that case, it needs to be chosen indepen- 
dently of /i 7 . 

Their theorem gives a concrete setup in which the KPZ formula holds and enabled 
Duplantier and Sheffield to state the following striking conjecture. 

CONJECTURE 1.12 (Duplantier, Sheffield |DSllj ). — Ifm n E M£ are sampled accord- 
ing to a statistical physics model in the 7 -universality class, then the pushforward mea- 
sures fi n of any "natural" embedding m„ C § 2 of m„ weakly converge asn->oo towards 
a random measure which is closely related to the Liouville measure /z 7 = e 7h , where h 



is an instance of the Gaussian Free Field on the sphere S> 2 ( n ) , See remark 2.20 for a 
definition of the GFF on the sphere § 2 . 

The reason why the limiting measure is not given exactly by the Liouville measure 
e lh stands from the fact that the measures \x n are renormalized to have measure one, 
while e lh has a random a.s. finite total mass. The actual limit is not given simply 
by conditioning e lh to have measure one, nor by renormalizing by J e' yh , it is slightly 



11. If one keeps track of the root, asymptotically it will be distributed according to /i 7 = e lh . 
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more subtle than that. See section 6 in [ShelO] for a precise conjecture on the limiting 
measure. 

The rest of this survey is divided as follows. In section |2j we will give a short 
introduction to Gaussian Free Field. Then, along section |3j we will give a meaning to 
the Liouville measures /i 7 = e lh . In [DSllJ, the Liouville measures are defined for all 
7 G [0, 2), we will give here a simplified proof which holds only for 7 G [0, \/2). Finally, 
in section [4], we will sketch the main ideas behind the proof of the main theorem. 

Acknowledgments: 

I wish to thank Vincent Beffara, Nathanael Berestycki, Cedric Bernardin, Nicolas 
Curien, Bertrand Duplantier, Gregory Miermont, Remi Rhodes, Mikael De La Salle, 
Scott Sheffield, Vincent Vargas and Wendelin Werner for very useful discussions and 
Bertrand Duplantier, Nam-Gyu Kang, Maxim Krikun and Scott Sheffield for some of 
the nice pictures that appear here. 

2. THE GAUSSIAN FREE FIELD (GFF) 

We will give here a short and by no means self-contained introduction to the Gaussian 
Free Field. We refer to |Dub09j IShe07j for complete references on this topic. We will 
try to give a certain flavor of what the GFF is and at the same time, we will introduce 
some of its key properties which will be needed later. 

2.1. Discrete Gaussian Free Field (DGFF) in the square 

Let us start by a discrete version, known as the Discrete Gaussian Free Field. To 
simplify, we will consider the case of the square domain [0, l] 2 . 

Definition 2.1. — The DGFF in Ajy := ^Z 2 n [0, l] 2 with Dirichlet boundary con- 
ditions is a probability measure on functions : An — > K such that h\g\ N = (where 
dA N : = ^Z 2 n d[0, l] 2 ) and with density 

dP[h N ] := - exp(-^ ^2(h N (x) - h N (y)) 2 ) Y[dh N (x) , 

where the sum is over nearest neighbor pairs x ~ y G A^r, and Z is a renormalizing 
constant to make it a probability measure. See figure [7] for a sample of a DGFF. 

Due to the Dirichlet boundary conditions, the density is also equal to Z~ l exp (l/2(tiN, Ah 
where A is the discrete Laplacian. Written this way, one sees that DGFF is a Gaussian 
random surface with covariance matrix given by A -1 . It is a standard fact that A -1 
is given by the matrix [G n{x , y)} x , y &A N , where Gtf(x,y) is the Green's function for the 
Random Walk in killed on dA^- (A good reference for discrete Green's functions is 
for example the book |LL10] ). As such, one may give the following equivalent definition 
of DGFF: 
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Figure 7. A sample of a DGFF on A 2 i. Picture by S. Sheffield. 



Definition 2.2. — The DGFF on An with D.b.c is the centered Gaussian process 
h 7v indexed by the points x G An and with covariance structure given by 

Cov[h N (x),h N (y)] = E[h N (x)h N (y)] :=G N (x,y). 

Remark 2. 3. — Just to give an idea of the amount of fluctuations, "in the bulk", say at 
x N ■= (N/2, N/2), one has 

(9) Var [h n (x N )] = G N (x N , x N ) x log . 

Remark 2.4- — If D ^ C is any domain of the plane, one can define similarly a DGFF 
/i/v with D.b.c. in the domain D by using the discretization Dn : = ^Z 2 fl D (and with 
boundary defined as dD^ '■= j^^ 2 H D c ). 

It turns out that as A oo, converges (in a certain sense to be precised later) 
towards a conformally invariant object called Gaussian Free Field. 



2.2. A first attempt at defining the Gaussian Free Field 



If we are following Definition 2.2, it is tempting to define the continuous limit of h 



N 



in the domain D as the centered Gaussian process indexed by the points x G D with 
the following covariance structure: 



(10) 



CovjTt(x), h(y)] := Go(x,y) for all x, y G D , 
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where Gd{x, y) is the Green's function of the domain D. See the later subsection 2.4 for 
a definition (and more) on continuous Green's functions. Unfortunately, this definition 
is "ill-posed" since, as suggested by Remark |2.3[ one would have for any x G D 



Var[/i(x)] := Gd(x, x) = oo . 

One way to overcome this would be to regularize h using smooth functions, and this 
is what we will do when we will introduce the e-regularization h e (z). Before, let us 



follow a different approach to define GFF inspired by our initial definition 2A_ 
2.3. GFF as a Gaussian process indexed by Sobolev space H 1 



Following Definition 2.1 it is natural to look for a probability measure on functions 
h : D — > R satisfying h\Qz> = whose intensity, informally would be given by 

(11) F[h] ocexp (- I [ \\Vh\ 

In order to find a well-defined object corresponding to this informal definition, it is 
useful to introduce the following Hilbert space: 

Definition 2.5 (The sobolev space — Let C = Cd be the set of smooth functions 
/ with compact support in D. We define the space H 1 (= 1-L\(D)) as the closure of C 
for the norm 

V/ 



2 _ ' / NT- /'II 2 
Z7T 



D 



(12) 



|2 



H 1 is a separable Hilbert space for the scalar product 

We are thus trying to define a random function h E H 1 such that 

P[/i] oc exp ^ — — \\h\ 

To gain some intuition, if the Hilbert space "H 1 happened to be finite dimensional, it 
would be some (IR fc , || • H2), with k > 1. In that case, if e±, . . . ,ef- is any orthonormal 
basis of M. k , then a random variable x6R fc with intensity P[x] oc e -1 / 2 " 1 " 2 can always 
be written 

k 

X ^ G>i e^ , 
1=1 

where (oj) are independent Gaussian random variables. 

Since 1-L 1 is a separable Hilbert space, let (e n )„>i be an orthonormal basis for T-L 1 . 
By analogy with the finite dimensional case, we want to define the Gaussian Free Field 
h as 

(12) h:=^2a n e n , 

n>l 



12. We do not pay attention to the constant in the exponential through this informal discussion. 
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where (a n ) are independent Gaussian variables ~ A/"(0, 1). The difference with the 
finite-dimensional case is that for any k > 1, the above Gaussian random variable x 
was almost surely in IR fc , while in our present case, it can be shown that almost surely, 
the above formal series (12) does not converge in TL l . 

In fact, it does not even converge in L 2 (D), and as such the Gaussian Free Field h 
will not be defined as a proper function, but instead as a generalized function in T>' (i.e. 
a Schwartz distribution). More precisely, it can be shown (see |Dub09j ) that almost 
surely, the above sum h converges in the space Hr x 



(13) 



Definition 2. 6. — From now on, the Gaussian Free Field with Dirichlet b.c. in a domain 
D will be defined as the random distribution 



h := a n e n a.s in H 1 



n>l 

where (e n ) n is an orthonormal basis of % 1 {D). (The Dirichlet b.c. is hidden in the fact 
that any / G H 1 satisfies f\go = 0). 

Example 2. 7. — In the case where D is the square [0, l] 2 , one can write down an explicit 
basis for W}([0, l] 2 ), namely for all j, k G N*, let 

(13) e jth (x, y) := = 2v / 2vr sm(jnx) sm(kny) . 

y 3 + k 2 

It is not hard to check that (ej t ).)j,ke^* i s indeed an orthonormal basis for (TL 1 , \\ ■ ||v)- 
A Gaussian Free Field in the square [0, l] 2 with zero boundary conditions can thus be 
written as 

h = 52 /-o' k ,o sm(jirx) sm(kTry) , 

'2 + k 



where (aj^)j,k are independent Gaussian variables of variance one and where the con- 
vergence for this series holds in the space T-L~ l . 

Definition 2.8. — If / = ^ n>1 ct n e n e "H 1 , then with a slight abuse of notation we will 
denote by (h, /) v the following quantity 



It is straightforward to check that for any / G Ti 1 , (h,f)^ is a Gaussian variable. 
More precisely, the following proposition follows easily from the definition of (h, /) v . 

PROPOSITION 2.9. — Let h be a Gaussian Free Field in D. Then the process 
((^j /)v)/e'H 1 ^ s a cen ^ ere d Gaussian process indexed by the set H 1 and with covariance 
structure 

cov[(/i,/) v , (M) v ] = (/>^)v f or an y f>9 e H- 1 . 



13. In fact, it turns out that a.s. h G H e for any e > 
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Remark 2.10. — In fact, this proposition can serve as another way to introduce the 
Gaussian Free Field. This is for example the point of view in [Dub09, She07j, where 
they introduce GFF as this Gaussian process indexed by H 1 . 



This approach thus gives a good generalization of definition |2.1 to the continuous 



setting. In fact it can be proved that the DGFF hjy weakly converges in the sense 
of distributions (for example in T-i^ 1 ) towards X-p. h, where \j2 is a lattice-dependent 
constant. 

The following proposition relates the Gaussian Free Field we have just defined with 
definition O 

PROPOSITION 2.11. - - Let h be a Gaussian Free Field in D. For any p G C D (the 
smooth functions with compact support in D ), we will denote by (h, p) the distribution 
h tested against the smooth function p. Then the process {(h, p)) c is a centered 
Gaussian process indexed by Cd with covariance matrix 



Gov[{h,p),(h,f/)] = // p(x)p'(y)G D (x,y)dxdy , 



DxD 



where Go is the Green's function of domain D. 
2.4. The Green's function Go of a domain D 

Definition 2.12 (Green's function in the domain D). — The Green's function of a do- 
main D ^ C will be denoted by the function Go : D x D — )■ R + . For any x G D, define 
the function G x (y) := Go{x,y). With such notations, the Green's function Go(x,y) is 
characterized by the following properties 

(a) AG X (-) = on D \ {x}, namely it is harmonic in D \ {x}. 

(b) G x (y) ->• as y -»■ 3D. 

(c) G x (y) ~ log as y -»■ x. 

By removing the logarithmic singularity, one can rewrite G x (y) in the following way 
G x (y) = log + G x (y), where G x (y) is the harmonic extension to D of the function 
— log on the boundary dD. 

Here are some well-known properties of Green's functions that we will use. 

PROPOSITION 2.13 (Properties of the Green's function). — 

(a) Conformal invariance: if <fi : D — >■ D' is a conformal map, then for any x,y G D, 

G D (x, y) = Gjy{4>{x), <t>{y)) ■ 

This follows easily from the definition of Green's function. 

(b) Note that G B (0,y) = log^. 

(c) Grj(x,y) = Go(y,x) (this can be seen for example using (a) with (b)). 

(d) For any x G D, G x eU 1 . 

(e) In the sense of distributions, ^A[G X (-)} = 5 X , the Dirac point mass at x. 
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(f) The above harmonic correction G x (y) satisfies 

G x (x)=\ogC(x,D), 

where C(x,D) is the conformal radius of D viewed from x. If (p is a conformal 
map D — > D with <f>(x) = 0, then C(x,D) is simply defined as This 
property (/) can be easily checked using (a) and (b). 



Let us now explain how one can recover (at least formally) Proposition 2.11 from 
the above property (e). For any / G Cd, one can make sense of ^-(Vh, V/) where V/i 
is understood in the sense of Schwartz distributions. It is not hard to check that this 
quantity is exactly (h,f)y. Now, since / has compact support, integration by parts 
implies 



Using this identity with p := [—A]/ G Cd, we find 



±(Vh,Vf) = ^(h,[-A]f). 



{h,p) = (h,[-2irA-']p) 



This implies Proposition 2.11 



since 



Cov[<M, (h,P')} = Cov[</», [-2 7 rA- 1 ]p>v> (K [-2^^%] 



([-27TA-V, [-27rA-V>, 
1 



(V[-27rA- 1 ]p,V[-27rA- 1 ]p / ) 

p(x)p(y)G D (x,y) dxdy 



2vr 

(p, [-2ttA-V> 



DxD 



where in the last equality, we used property (e). 

Finally, let us mention that using property (b), it is not hard to extract the following 
striking property for Gaussian Free Field. 



PROPOSITION 2.14 (conformal invariance). — Let <p : D' — )■ D be a conformal map. If 
h is GFF in D, then h! := h o <f> is a Gaussian Free Field in D' . 



2.5. The e-regularized GFF h e 

The purpose of this subection is to regularize the Gaussian Free Field h in order to 
obtain a smooth function h e . For this, we will rely on the following e-regularization of 
the Green's function. For any e > and any point x G D, let 

G X M :=log— p ; + G x (y). 

e V \x — y\ 

This regularization has the following important property: 
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PROPOSITION 2.15. - - For any x G D, and any e > 0, G x e G H 1 . Furthermore, in the 
sense of distributions, one has the following identity 

(14) ^A[G?(.)] = iv, 

where v x ^ denotes the uniform measure on the circle of radius e around x, dB e {x). 

Remark 2.16. — In fact, with the above definition of G x , the proposition as stated is 
not correct when x is close to the boundary (d(x,dD) < e). To overcome this issue 
while keeping the same statement for the proposition, the definition of G x has to be 
modified accordingly near the boundary dD. To keep things simple, we choose in this 
paper to neglect these effects. We refer to |DSllj where this technicality is handled 
properly. 

This regularized Green's function enables us to introduce h e (z) the GFF evaluated 
against u XjC , the uniform measure on dB e (x). Informally, it corresponds to h e (z) : = 
(h, v x ,e)- L e t us define it as follows 

Definition 2.17. — If h is a sample of a GFF in D, then for any z G D, let 

h € (z):= (h,G*) v , 
which is well defined since G x e e % x . 

In fact, it corresponds exactly to our informal definition thanks to the following 
computation: since G* G H , one has 

h e (z) = (h,G*) v 

= ^-(h,[-A]G:) 

= (h, V x ,e) 

2.6. A Brownian motion out of the GFF 

Along this subsection, we will identify a very useful Brownian motion "within" the 
Gaussian Free Field. We start with the following lemma: 



LEMMA 2.18. — For any z G D and any e > < - 14 - ) , one has that 



(15) Vai[h e {z)] = log^ +logC{z,D) 

where C(z,D) is the conformal radius of D viewed from z. 



14. To be self-contained here, one should assume here that d(z,dD) > e 
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Proof: 

Var[/i e (z)] =Var[(/i,G*> v ] 
= (G z ,G e z ) v 

= log- + J G z (x)dv z>e (dx) 

= log* +logC(z, J D)- □ 
The following proposition will be of crucial importance in the remaining of this text: 

PROPOSITION 2.19. - - Let h be a GFF with zero-boundary conditions in some domain 
D. For any point z e D, let £q := mi{t > : B e -t(z) C D} and let 

Y t (z) := h e -t(z) , 

be the stochastic process defined for any t > (Recall h € (z) denotes the above regular- 
ization). 

Then with such notations, the stochastic process 

Bt(z) :=Y t * +t -Y t >, 
is a standard Brownian motion. 

Proof: The family of random variables {B t (z)} t > is clearly a Gaussian process. 
Therefore, it only remains to check that for any < s < t, Cov[B s (z),B t (z)] = s(— 
s A t). Let r := e~*°, r\ := e~^° +s ^ and r 2 := e~^° +t \ Let us first compute 

Cov [h ri (z) , h r2 (z)] — E [(h, G z ri ) v (h, G z 2 > v ] 

' G r2 ) y 

= \og^^ + \ogC(z,D) 
= t z + s + \ogC(z,D). 

One can compute in the same way Cov[h ro (z), h ri (z)\ and Cov[h ro (z),h r2 (z)\. This 
gives us 

Gov[Bs(z),Bt(z)] = Cov[h ri (z) - h ro (z), h r2 (z) - h ro (z)] 

= t Q + S + t Q 2 * £q 

= s □ 

Let us conclude this section on the Gaussian Free Field by the following remark. 
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Remark 2.20 (Gaussian Free Field on § 2 ). - - We have just defined the Gaussian Free 
Field on a domain D with Dirichlet Boundary conditions. In the same fashion, one can 
define a Gaussian Free Field on the sphere S 2 (this is needed for example if one wants 



to make sense of Conjecture 1.12). In this case, the Green's function is given by 



:= log [cotan-] for all x, y G § 2 , 
where 9 denotes the angle between x and y. If one wants to define the Gaussian Free 



Field on S 2 as a Gaussian process similarly as in definition 2.9 , the construction can be 
done in the same fashion except that in this case, the natural Hilbert space to consider 
would be the closure for the norm || ■ ||y of the space C^(S> 2 ) of smooth functions 
<p : § 2 — > M. with L 2 (ft(x)dx = (the integral here is with respect to the area measure 
on§ 2 ). 



3. THE LIOUVILLE MEASURES e 7?i 



The purpose of this section is to make sense of these Liouville measures e 7/l , which are 
crucial in the main Theorem 1.10 The approach followed in |DSllj is to discretize e lh 
into e lhf - (where h e is the e-regularization of the GFF h we have introduced in subsection 



2.5) and to then let e — » 0. As it will become clear below, without renormalization, e 



7 he 



would diverge in the space of measures. The natural discretization will be the following 
one: 



Definition 3.1. — For any domain D, any 7 > and any e > 0, let /x e be the measure 
absolutely continuous with respect to Lebesgue measure C and such that 

(16) dfi t {z) := e^e fh ^ z) dz. 

Duplantier and Sheffield prove the following proposition in |DS11) : 

PROPOSITION 3.2. — J/7 G [0,2), then for any domain D, almost surely as e \ 
along powers of two, the measures fi t weakly converge inside D towards a non-degenerate 
random measure /x 7 which we will call the Liouville measure of parameter 7. The 
Liouville measure /x 7 is measurable with respect to the Gaussian Free Field h and we 
will denote it sometimes by /x 7 = e 7 h . 

Remark 3.3. — 

(i) If 7 G (0,2), it can be shown that the Liouville measure /i 7 is a.s. singular with 
respect to Lebesgue measure. 

(ii) If 7 > 2, in some sense things become "singular". See for example the work 
[BJRV12] which studies this case. 
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We give here a new proof of this proposition which holds only for the regime 7 G 
[0, \/2) and furthermore our convergence result will hold only along a certain subse- 
quence \ that we will not make explicit (in [DSllJ, it is also along a particular 
subsequence, but they show that e& = 2~ k is enough). We believe this proof is interest- 
ing in its own since it is slightly different as the one carried in jDSllj . yet it cannot be 
extended to the range 7 G [y/2,2). See the proof in jDSTT] which gives the full range 
[0,2). 

Proof in the case 7 G [0, a/2): 



To simplify, we will restrict ourselves to the case where D is a bounded domain. We 
wish to prove the following proposition: 

PROPOSITION 3.4. — J/7 G [0, y/2), then for any continuous function : D — > R, the 
sequence of random variables 



{/^0)} £ >o = \ / <t>(z)^e(dz) 

'D J e >0 



is a Cauchy sequence in L 



2 



3.2 



in the regime 7 G [0, \/2). 
t is well-known that this 



Let us first see why this proposition implies Proposition 
Let QJt(Z)) be the space of finite positive measures on D. 
space equipped with the weak* topology ( called "weak convergence of measures" in 
Probability theory) is a complete, metrizable, separable space. Here is an example of a 
metric on VJt(D) which induces the weak* topology. Let (4>j)j>i be a countable basis of 
the separable space (C°(D), \\ ■ ||oo) of continuous functions on D such that ||0j||oo < 1 
for all j. Then 

defines a metric on VJl(D) for the weak* topology. 



Using Proposition 3.2, one can find a subsequence (ek)k>i such that uniformly for all 
j G {1, ...,*}, 

E[(/i ev (^) - ^ fc „(0,)) 2 ] < 2- 3fc for all k" > k> > k . 

By Markov's inequality, this implies P [|/4 fc , (<pj) — fi ek „ (<pj) \ > 2 _fe ] < 1~ k for all k > 1 
and j < k. Using Borel-Cantelli lemma, it is an easy exercise to show that this in turn 
implies that fi tk is a.s. a Cauchy sequence in (SOT(Z)), d). Since the later space is 
complete, we thus obtain an almost sure limit /1 = /i 7 G VJt(D). Furthermore, since 
each random measure fi ek is clearly measurable with respect to the Gaussian Free Field 
h, we obtain that /i 7 = lim H tk {h) is itself measurable with respect to h as a limit in 
(5DT(-D),gQ of h- measurable measures. 



Proof of Proposition 3.2 Let us start with the simpler lemma 
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LEMMA 3.5. — If 7 £ [0, \/2), then for any continuous function <fi : D — )■ R ; we /iave 
E[^(0) 2 ] — ► // <f>(x)<t>(y)[C(x,D)C(y,D)Y a/ \^ G »^dxdy. 

JJdxD 

Proof of the lemma: 

E[/i e (0) 2 ] = /Y 0(x)<P(y) ^ E[e^ {x)+ ^ y) ]dxdy 

JjDxD 

(17) = // ^(x)(f)(y)e^ e^ Yar \- he{x)+he{y) hxdy , 

J JDxD 

where we used the fact that if X ~ A/"(0, cr 2 ), then its Laplace transform is given by 

E[e^ x ] =e*£. 
Now, 



V&r[h e (x) + h e (y)] = Var[h e (x)] +Vai[h e (y)] + 2 Cov [h e (x), h t 
= Wsr[h e (x)] + Var[h e (y)] + 2 (G x , v e , y ) . 
If I a; — y\ > e, we find exactly 



(15) 



Var [/i e (ar) + /i e (y)] = 2 log - + log C7(ar, £>) + log C(y, D) + 2 G D (x, y) . 

Let H(x,y) := log[C(a;, D)C(y, D)] + 2Gi)(x,y), then if |a; — y| < e, one finds instead 
the inequality 

Var [h e (x) + h e (y)] < 2 log - + H(x, y) . 



Plugging these into (17) gives us 



E^(0) 2 ] = [[ ( p(x)4>(y)[C(x,D)C(y,D)f / \^ G ^dxdy 

J J\x-y\>e 

' H{x ' y) dxdy) 



1 
e 2 



+ 0(\\<f>\L 

'\x-y\<e 

To conclude the proof of the lemma, one needs to show that the second term goes to 
zero as e — > while the first one remains bounded. The key contribution in both cases 
is what happens when x ~ y. In that case, we know that Go(x,y) ~ log ,_, . This 
implies that when x ~ y, the term e' y2GD ^ x,y ^ behaves like 



1 2 G D {x,y) 



x-y 



7 2 +o(l) 



where o(l) — > as x — > y. Using the fact that if a < 2, then JJ DxD 
an easy exercise to conclude the proof of the Lemma. 



x-y 



< 00, it is 
□ 



15. in fact we also need to assume here that d(x, 3D) A d(y, 3D) > e, but we neglect these boundary 
issues here which are easy to be taken care of 
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Now, let us prove Proposition 3.4, i.e. that {/x e (0)} e> o is a Cauchy sequence for any 
continuous on D. For this, let us estimate for < rj < e: 

E[(/i e (0) -M0)) 2 ] =E[/i e (0) 2 ] +E[^(0) 2 ] 

-2e^rj^ [[ 4>(x)4>(y)e^ YaT l- h ^ x)+h ^ y) \dxdy. 



J JDxD 

Similarly as in the proof of the Lemma, we find that 

Var [h e (x) + h v (y)] < log - + log - + log(C(x, D)C(y, D)) + 2 G D (x, y) , 

with equality if and only if \x — y\ > e V 77 (and <i(x, <9.D) A c?(y, <9.D) > e as well). In 
particular, in the same fashion as above, if 7 is chosen so that 7 < a/2, this implies that 



e 2 77" 



DxD 



0(ar)0(y) [C(x, D)C(y, D)] 7 ^ 2 ^ 2 ^^ dxdy . 



D 



0<77<e->0 JJ DxD 

which thus implies that /x e (</>) is indeed a Cauchy sequence in L 2 . □ 

Remark 3.6. — Another natural approach would be to consider e 1— > yU e (0) as a stochas- 
tic process in e \ 0. Since it can be written as 

7^(0) = f <j)(x)e^ {x)+ ^ lose dx, 

(j){x)C{x,D)-^e lK{x) -^ Y&r \- K{x) ~\ , 

D 

and since for each fixed 1, gTeM-^-Var^O)] - g & p OS j^j ve martingale, one might 
be tempted to prove the a.s. convergence of /4 when > 0, by showing that 
it is a positive martingale (furthermore the above L 2 bounds when 7 < \/2 would 
imply its uniform integrability). Unfortunately, this is not the case. In [DSllJ, the 
authors manage nevertheless to rely on such an approach by looking at a different way 
to regularize h, namely h n := ^2 i<n {h, /«) v /j, where {/j} is some orthonormal basis 
of ft 1 . In that case l^n{4>) (defined accordingly, see |DS11] ) is in this positive 
martingale. 

Remark 3.7. — Finally, one should point out that such measures had already been 
constructed within the theory of Gaussian Multiplicative Chaos initially developed 
by Kahane in |Kah85] . See [RV10J for a more general construction of such measures. 



4. IDEAS BEHIND THE PROOF OF THE MAIN THEOREM 



In this section, we wish to explain where the KPZ formula comes from by giving some 
of the ideas behind the proof of Theorem 1.10 
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4.1. Setup 

Let us fix some (deterministic) K C [0, l] 2 and some parameter 7 G [0, 2). Assuming 
that the limit exists in equation ([7]), our goal is to express the quantum scaling exponent 
A = A (A) as a function of the Euclidean scaling exponent x = x(K). Recall that 

A = A(K) : = lim — ±± , 

s^o log 

where z is "sampled" according to the a.s. finite measure /i = e lh and E averages over 
the random measure [L. In this setting, we first sample /i 7 = e lh ~ E and conditioned 
on the measure /i 7 , we sample z ~ /i 7 . For the proof of the main theorem, it will be 
useful to invert this procedure, i.e. to first sample z (according to the correct marginal 
measure) and then to sample /i 7 conditioned on the value of z. The coming subsection 
introduces the right tool for this. 

4.2. Rooted Liouville measures 

As we have just explained, the definition of the exponent A involves the coupling 
(fj, y , z), where z is sampled according to the first coordinate /x 7 (which in general is not 
a probability measure). The law of this coupling can be written dh x dfj,^(z). In order 
to invert this sampling procedure, we will introduce the following probability measure 
which can be viewed as an e-regularization of the above coupling. 

Definition 4-1 (Rooted Liouville measure). — For any 7 G [0, 2) and any e > 0, let 

6 e := ^e^ {z) dhdz , 

be the probability measure on 1-L~ l x [0, l] 2 , where Z e is a renormalizing constant chosen 
so that e is a probability measure. 

This regularized coupling enables us to make sense of the reversed sampling proce- 
dure. First of all, we need to compute the marginal distribution of e on z. It is simply 
given by p e (z) = Z~ 1 ¥ l h[e" /h ^ z ^~\ which is explicit since using Lemma 2.18 one has that 
it is proportional to C(z, U) l2 l 2 when d(z,dD) > e. In particular, as e — > 0, the den- 
sity p e (z) converges towards a limiting p(z) oc C{z^D) 1 I 2 . Now, conditioned on z, the 
marginal on h is given by the Gaussian Free Field h weighted by e jhe ^ i.e. by j^ht^dh ' 
Note that this step would not make any sense without our e-regularization since the 
scalar quantity "e 7?t ( z )" is not defined at point z. It is a standard fact about Gaussian 
processes (Cameron-Martin theorem), that if h z is sampled according to j^lt^l)^ ? then 

h z = h + 7(7*. To see this at least heuristically, note that the process h + 7 G\ is a 
deterministic translation of the Gaussian process h and it is easy to check that in the 
finite dimensional case, if A is a standard Gaussian vector in R n , then for any fixed 
u G lR n , the law of X = X + u is the same as the law of A weighted by e( x ' u K In 
the case of the Gaussian Free Field, the scalar product is (•, -) v and we indeed have by 
definition that (h,jGf) v = •yh 6 (z). 
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In particular, the pair (z, h z ) ~ B e can also be sampled as follows: first sample z 
according to the above marginal distribution p e (whose density away from 3D is pro- 
portional to C(z,D)' y / 2 ), and then conditioned on z, let h z := h + 7G 2 , where h is 
an independent Gaussian Free Field in [0, l] 2 . This way, we see that the measures e 
converge towards a limiting measure as e — > 0, for which the GFF (h z ) conditioned 
on the first component z ~ p is sampled according to h z := h + 7 G 2 . 

Using this coupling we have that 

E[p[B s (z)nK ^ 0]] x e\B 5 (z) HKy£ 
(18) x f p(z)dz ¥ h [Bl =h+lGZ {z) n K ± 0] , 

</[0,l] 2 

where Bf lZ=h+ ^ GZ (z) denotes the quantum ball of quantum area 6 around z in the sense 
of Definition 1.8| but with a field sampled according to h z := h + 7 G z . 



4.3. An estimate about quantum balls around the root 



For any z G [0, l] 2 , in the same fashion as in Proposition 3.2, it can be shown that if 
7 < 2 and if h z := h + 'y G z , then one can make sense of the Liouville measure rooted 
at z, p z := e lh \ 

In [DSllJ, the following property is shown. 

PROPOSITION 4.2. — Let z e [0, l] 2 and let h z ■= h + ^G z , then almost surely as the 
Euclidean radius r \ ; 

where Q = Q 7 := 2/7 + 7/2 > 2 and where h z := h r + "iG z r . (c = c 7 is some explicit 
constant). See subsection J^.l in |DS11| . 

We will not prove this proposition, but instead we will convince ourselves through 
the computation of the expectation of /i 2 (i? r (z)) that one can indeed expect such a 
behavior: 

Lemma 4.3. — 

E[n*(B r (z))] ~ cr 7Q E[e 7 ^ {2) ] , 
as r — > for a certain constant c = c 7 . 



Proof: First of all, using Lemma 2.18, one has (if d(z,dD) > r): 

= C(z,D)£r-£ C(z,Dy 2 r^ 2 
= C{z,D)^r-^ 2 ' 2 . 
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In order to compute E[//*(S r (z))] , let us approximate h z = h + 7G 2 into h\ 
h e + jG z : 

E[^(B r («))]=Jme^ f ^[e^ {x)W G * {x) ] dx 

J B r 

= / e^" 6z(x) e^^^dx using Lemma |2. 18 



e \G z (x) e 7 2 \og\x- z^+^G^x) 



f exp [ 3 ^G z (z) + o(l)] r ^—,dx 

JB r 2 SB - 2| 7 



as r — > 0, since x 1— >■ (5* (a;) is continuous. Therefore as r — >■ 0: 



E[^(s r (z))]~c(z,£>) 3 ^ r -^—^dx 

JB r \ X ~ Z \ 



u=0 



Remark 4-4- — This first moment computation indeed provides some supporting evi- 
dence for Proposition 4.2| Yet, such a comparison of first moments is not so natural 



after all since, as it was pointed out to us by Nicolas Curien, the expected quantum 
area of B r (z) diverges as r — > when 7 E (y/2,2) (in the above displayed equation, 
2 — 7 2 < when 7 > y/2). This counter-intuitive phenomenon is due to the fact that 
for any 7 > 0, the main contribution in E[e 7 ' l ^ z - ) ] does not come from typical prop- 
erties of h e (z) but follows instead from large deviations events for h e (z). This is why 
first moments computations are not suitable for studying "typical" behavior as one is 



interested in in the statement of Proposition 4.2 See subsection 4.1 in |DS11] for a 



proof of this "law of large numbers" type of behavior. 



Recall that the quantum ball B 5 {r) for a field h z is defined as B T (z) with r 



sup{r > 0, y!t{B r {z)) < 5}. The content of Proposition 



4.2 



tells us that B iz) should 



be very well approximated by the ball B (r) defined as Bf(z) with 



r := sup [r > s.t. cr 



<S}. 
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Plugging this into (18) gives us 



E[fi[B 5 {z)nK ^0]] 



[o,i] 2 



[o,i] 2 



[o,i] 2 



[o,i] 2 



p(z)dz¥ h [B i 



h z =h+~fG z 



n K ^ 0] 



Eft 



lB f (z)nii"^0 



As we will see below, it is not difficult to show that the law of the random radius 



f = f(z) at point z depends very little on the point z G [0, l] 2 ( 16 - ) . In particular, if P 



denotes this common law for r, we have that 



E[fx[B s (z)nK^0\]tt eflP(f) / dzl Br{z)nK ^ 

JR+ ^[0,l] 2 

« / dP(f) f 2x{K) by definition of x — x(K) 
Jr + 

Therefore it only remains to understand the law of f (in some sense uniformly in the 
root z G [0,1] 2 ). This will be done by identifying a drifted Brownian motion within 

h* :=h + "fG z . 



4.4. Reduction to a large deviation question on Brownian motion 

Let us fix some z G (0, l) 2 and let to > so that B e -t (z) C (0, l) 2 . Then similarly 



as in Proposition 2.19 if 



Wt(z) 



hi 



— h z 
t n e - t 

K-*-*o ( z ) ~ K-*o (z) + 7 [G z e -t-t 0) 
h e -*-t (z) - h e -t (z) + 7 1 , 



G z e - t0 (z)] 



16. we are neglecting boundary issues here 
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then (W t (z)) t >o is a Brownian motion with drift 7. I.e W t — B t +~/t. r can be defined 
using this Brownian motion: indeed recall 

f := sup {r > s.t. C r 7Q e 7 ^ (2) < 5} 



exp — 
exp — 

exp — 
exp — 



inf U >0 s.t. c(e~ t V Q e' yh ^ {z) < 5 



}] 



inf 1/ > s.t. ( e -*)^e 7Bt+72 * < 5}] 

log 5 



inf a > s.t. B t + (7-Q)t< 



7 



7 



inf <U > s.t. B t + a^t> 



where B t := (—Bt) is a standard Brownian motion and where a 7 := Q 7 — 7 = 2/7 — 7/2, 
which is positive when 7 < 2. Let T = Ts be the stopping time for the drifted Brownian 
motion B t + a 1 t stopped the first time it reaches level 7 _1 log|. Since f w e~ Ta , 
summarizing the above discussion, we obtain that A = A(K) should be given by 

logE[e" 2a:T5 ] 



A = lim 



log 5 



It remains to compute the quantity E[e 2:cT ' 5 ]. This is a classical computation for 
Brownian motion and it works as follows. Consider for any /3 > the process: 

6 2 

t ^ exp (pB t - yt) , 

which is a martingale. Using the optional stopping theorem for the stopping time Tg, 
we get for any (3 > 0, 

E [exp ( - a 7 T s + /3/ 7 log ^ - /3 2 T s /2) ] = 1 , 

which in turn gives 

E[e" 2 ^] =5^\ 

if ft — /3 7 is chosen so that 2x = /3a 7 + /3 2 /2. Since A = A(K) is given ultimately by 
/3 7 /7, this indeed gives us a quadratic relation between x(K) and A (If). One can check 
that this quadratic relation is the KPZ formula (pj). □ 



4.5. Other proofs of a KPZ formula in the literature 

Finally, let us mention that after Duplantier and Sheffield announced their result, 
other proofs of KPZ formulas have been proved in slightly different settings: 

- Benjamini and Schramm obtained in |BS09] a simple and enlightening proof of a 
KPZ formula for multiplicative dyadic cascades in dimension one. The advantage of 
their proof is that it gives a quadratic relation between actual Hausdorff dimensions 
as opposed to "expected box-counting dimensions" in [DSllJ Unfortunately, their 
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argument is inherently one-dimensional and if one would extend their argument to 
higher dimensions, it would no longer deal with proper Hausdorff dimensions. 
Rhodes and Vargas proved in [RV08J a KPZ formula in the general setting of Gaus- 
sian multiplicative Chaos. Their proof enables to deal with "stationary" measures 
(as opposed to [BS09J which relies on a discrete dyadic division). The difference 
with [DSllJ is that their KPZ formula holds for a different notion of dimension (or 
rather scaling exponent) as the one considered in |DS11] . In that sense their work 
is complementary to the work |DSllj . More precisely, in rough terms, if K C [0, l] 2 
and if \i denotes a measure (for example the Liouville measure), then their notion 
of scaling exponent is defined as 

A(tf):=infise(0,l],s.t. inf {V/^B^r,))*} = 

I coverings KGUB(xi,ri) ' ' 

where the balls B(xi, r») are Euclidean balls of radii rj. This notion is very different 
from the expected box-counting dimension considered in jDSllj . Note that in both 
works jDSllj and |RV08j. the notions of scaling exponents A, A still rely somewhat 
on the Euclidean metric. It seems one is still far from a "true" KPZ correspondence 
between Euclidean and quantum metrics. 
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